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Abstract 

The application of real space renormalisation group methods to quantum lattice models has 
become a topic of great interest following the development of the Density Matrix Renormalisation 
Group (DMRG) by White. This method has been used to find the ground and low-lying excited 
state energies and wave functions of quantum spin models in which the form of the ground state is 
not clear, for instance because the interactions are frustrated. It has also been applied to fermion 
problems where the tendency for localisation due to the strong Coulomb repulsion is opposed by 
the lowering of the kinetic energy which occurs as a result of electron transfer. The approach is 
particularly suitable for one, or quasi-one dimensional problems. 

The method involves truncating the Hilbert space in a systematic and optimised manner. 
Results for the ground state energy are thus variational bounds. The results for low-lying ener- 
gies and correlation functions for one dimensional systems have unprecedented accuracy and the 
method has become the method of choice for solving one dimensional quantum spin problems. 

We review the method and results obtained for the spin-1 chain with biquadratic exchange 
as well as the spin-1/2 model with competing nearest and next nearest neighbour exchange will 
be described. More recently, the DMRG has been applied to reformulate the coupling constant 
renormalisation group approach which is appropriate for the study of critical properties. This 
approach has been applied to the anisotropic spin-1/2 Heisenberg chain. 

Finally, we discuss recent work which has borne promising applications in two dimensions — the 
Ising model and the two dimensional Hubbard model. 



1 Introduction 



Recent times have seen the fabrication of an increasingly rich variety of materials whose low energy 
magnetic properties require intrinsically quantum mechanical as well as low dimensional lattice models 
for a realistic description. On the other hand, the development of sophisticated field theories and the 
surprising conjecture of Haldane has dramatically rekindled theoretical and mathematical interest 
in one dimensional quantum spin models. Also, the advent of high temperature superconductors has 
generated intensive interest in the magnetic properties of two dimensional quantum lattice models 
such as the Hubbard, t-J and Heisenberg models. 

It is well known that analytical solutions of these models are only available in rare and special cases 
and that mean field theories perform unreliably when applied to low-dimensional systems. Also, the 
presence, for example, of frustration in these systems renders the construction of suitable variational 
wave functions very difficult as the structure of the ground state is not at all clear ahead of time. 
As a result, perturbative attempts at systematically improving upon variational model states are 
often prone to slow or even spurious convergence. This problem of frustration has also for a long 
time plagued attempts at applying Monte Carlo methods to many of the abovementioned systems at 
low temperatures. Exact diagonalisation approaches may fail to reach suitably large lattices for two 
dimensional or coupled chain problems due to the exponential growth in the size of the Hilbert space. 
This is especially problematic for frustrated systems that possess incommensurate order in the ground 
state. 

There has therefore been a real window of opportunity open for the development of a portable, 
robust and systematically improvable numerical scheme for solving low dimensional quantum lattice 
models of magnetic systems. The advent of the renormalisation group (RG) and Wilson's highly 
successful numerical solution of the Kondo problem j|] in addition to promising applications to classical 
systems in two dimensions || suggested that the Real Space Renormalisation Group (RSRG) might 
well fill this role. 

Unfortunately, attempts at carrying out this programme produced results with disappointingly slow 
convergence The underlying causes for the poor performance of the RSRG were addressed in a 
seminal series of papers by White || where a powerful improvement on the RSRG, the Density Matrix 
Renormalisation Group (DMRG), was formulated. These papers were followed with a phenomenally 
successful study of the spin-1 chain where the Haldane conjecture, that spin chains of integral spin 
are gapped, was verified with unprecedented accuracy || . A flood of applications to problems in low 
dimensional quantum magnetism has since followed. 

In this article we review the DMRG and some its applications to problems in magnetism. In 
section ^| we briefly outline the method and the quantities of interest which can be calculated. In 
section [| we discuss applications of the DMRG to spin chains with frustration. In section [| we review 
attempts at using the DMRG to study phase transitions in these systems, either by direct calculation 
of order parameters or by deriving RG transformations in the space of coupling constants. In section |5| 
we discuss recent extensions of the DMRG to two dimensional problems and one dimensional models 
at non-zero temperature. 
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2 The DMRG method 



The DMRG is an iterative, truncated basis procedure whereby a large system (or superblock) is built 
up from a small number of lattice sites by adding a few sites at a time. As mentioned, the DRMG 
was introduced and formally developed in a series of papers j^] by White and co-workers at Irving. 
We refer the reader to these excellent papers for a detailed description of the method. Here we will 
attempt to present a simple exposition which captures the essential algorithmic steps involved. It will 
be useful to work with a concrete model; the one-dimensional, spin-1/2 XY model 

n = (s+sr+i + h.c.) (i) 

i 

where Sf 1 is the spin-1/2 raising (lowering) operator for lattice site i. 

Let us initially define a system block A which consists of only a single site and hence has M = 2 
states |1) = |t) and 2} — We record the matrix elements of the spin operator S A for this site, 
as well as the block Hamiltonian Ha, with respect to the system block basis {\n) : n = 1, . . . , M} i.e. 
S(n',n) = (n'| S% \n) and H(n',n) = (n'\H A \n) = 0. 

In the first iteration we form an environment block B which is precisely the same as the system 
block (a single site). Combining the system block, the environment block and two extra sites, we form 
a superblock A • •£>. A basis for the superblock is given by 

{\n s ) = Wa) <8> |si) ® |s 2 ) ® \n B ) :n A ,n B = l,---,M; si,s 2 =f, 1} (2) 

We write down the Hamiltonian 7is = Ti.A + T~Ca* + + H,b + "Hb for the superblock viz 

(n' s \n A \n s ) = 8 s > iSl 8 s > 2S2 5 n > BnB H(n' A ,n A ) (3) 
(n' s \H.A. \n s ) = S s ' 2 s 2 Sn' B n B [S(n' A ,n A ) (s[\ S, \sx) + h.c] (4) 

Low energy eigenstates of the 4-site Hamiltonian Tis may be found at this point. 

Next, we expand the system block by augmenting it with one of the extra sites. A basis for the 
augmented block A = A* is given by 

{ \n)) = \n A ) ® \s x ) :n A = l,...,M; a x =T, 1} (5) 

and clearly has M = 2M states. 

In order to again form a (larger, by two sites) superblock, we require the matrix elements of the 
augmented block Hamiltonian H A = TLa + Ha* as well as the spin operator corresponding to the 



added (end) site viz 

H(n',n) = ((n'\H A \n)) (6) 

= 5 s > iSl H{n' A ,n A ) + [S{n' A ,n A ) (s[\ S. \a x ) + h.c] (7) 

S(n',n) = «n'|S+|n» (8) 

= 5n >A (s'i\S:\ Sl ) (9) 
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The superblock construction (j|), (||) and (Q) can now be extended indefinitely where A i— ► A, 
H i — > H , S i — ► S , Al i — > M and |.)} i— > |.). However, even with the aid of sparse matrix diagonalisation 
methods, the size of the superblock basis, 0(M 2 ) will soon grow prohibitively large as M is doubled 
at each iteration. 

Instead, we truncate the basis of the augmented block so that the number of states used in the 
subsequent iteration does not exceed a fixed threshold, m. The truncated basis is chosen so as 
to contain the most important states for forming low energy eigenstates of the Tis- White's most 
important innovation was the use of a density matrix to determine these states. 

For instance, if we are interested in the ground state \4>g), then we form a reduced density matrix 
p for the augmented block as follows 

((n'\p\n))= <nsl0G>^G|ns) (10) 

S2 ,TIB 

with 

|»s) = \ n A) ® W) ® |s 2 ) ® \n B ) ee \n)) ® \s 2 ) ® \n B ) 

\n' s ) ee \n' A ) ® \s' x ) <g> |s 2 ) <8> \n B ) = \n')) ® \s 2 ) <8> \n B ) 

That is, the superblock ground state projection operator is traced over the environmental degrees of 
freedom uB. 

The density matrix eigenstates { \n))) : n = 1, . . . , M} are found by standard dense matrix diag- 
onalisation methods and the m most important states (corresponding to the m largest eigenvalues of 
p) are retained in forming the system block for the next iteration. Of course, the operators of interest 
must be represented in terms of these states viz 

H{n\n) ^ (((n'\n"))H{n" : n"')((n"'\n))) ^ H(n',n) (11) 

n"n"' 

where 1 < n, n' < min{M, m} i— ► M. 

The density matrix eigenvalues sum to unity and the truncation error, defined as the sum of the 
density matrix eigenvalues corresponding to discarded eigenvectors gives a qualitative indication as to 
the accuracy of the calculation as well as providing a framework for extrapolation to the m = oo limit 
H. As a result of basis truncation, the cpu time is linear in the superblock size and so one can easily 
proceed to the thermodynamic limit for intrinsic properties such as the ground state energy density 
(see jjj for a formal discussion of the DMRG fixed point and the thermodynamic limit of the DMRG). 

As mentioned, the accuracy of results obtained in this way was unprecedented [0], the accuracy 
of the ground state energy density eo for the spin-1 chain being ultimately limited by the precision of 
machine arithmetic viz eo = —1.401484038971(4). Similarly exquisite accuracy persists when targeting 
low-lying excitation energies in various symmetry sectors Q e.g. E\ — Eq = 0.41050(2). 

The method described above is the infinite lattice algorithm. In the finite lattice algorithm, the 
infinite lattice algorithm is initially applied so as to obtain a superblock of a desired size. System 
blocks are retained at all iterations and are used as environment blocks for the purpose of recursively 
deriving progressively superior system blocks within the superblock whose size remains fixed p|. 
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Repeated sweeps across the lattice result in substantial improvements in accuracy, especially for the 
study of edge states and systems with impurities. 

Applications to spin chains with high spin [B, ||, dimerisation and/or frustration |tTof— p^|, to- 
gether with extensions to coupled spin chains pi , models with itinerant fermions |L7]], |Q, Kondo 
systems [^9|-j2^], as well as coupled fermion chains [[23| have followed, even in cases of intermediate 
doping. Formulations for systems with single |24| as well as randomly distributed p5f impurities and 
disorder [^6| have also been forthcoming. 

By retaining matrix elements of operators acting on sites throughout the blocks, it is possible 
to calculate static correlation functions. In || and |^7j it is shown that highly accurate correlation 
functions can be calculated over many lattice spacings. Highly accurate studies of the structure factor 
and string order parameter (topological long range order) ||, as well as edge states Q in Haldane 
phase systems have also been performed. Dynamical correlation functions have also been calculated 
within the static |29| , continued fraction J30| and correction vector (3lJ approaches. Finally, the DMRG 
has been formulated for the solution of models of spin chains dynamically coupled to (dispersionless) 
phonons p2[ . 



3 Spin chains with frustration 

A spin model on a bipartite lattice with a nearest neighbour, antiferromagnetic Heisenberg interaction 
can have its tendency towards commensurate order in the ground state frustrated in the presence of 
other interactions. In this section we consider two such models in one dimension. In each case, J\ > 
denotes the strength of the Heisenberg interaction and J2 > the strength of the frustrating term. 

The Majumda-Gosh (MG) model |Q has the Heisenberg interaction frustrated by the presence of 
next-neighbour interactions viz 



H 



jyiSi.Si+i + j^ 10 ^ Si.Si+2 (12) 



where Si denotes the spin 5=1/2 operator for lattice site i. This model has some relevance to the 
interpretation of experiments on quasi one dimensional spin-1/2 chains J35|. For 5=1, the order 
can also be disturbed in the Lai-Sutherland (LS) |36| model through the presence of a biquadratic 
term 

H = E [ J i^+i + ^ i5) (5,.5 i+1 ) 2 ] (13) 



The models (|2|) |33j-|35|, |7| and (|13|) |3J, [|38|-|39[ have been the subject of intense study in 
the past. Exact diagonalisation methods together with exact solutions and rigorous results at special 
points have been combined to give a good picture of their critical properties. 

In the MG case, the pure Heisenberg point (J2 = 0) has been extensively studied |l0| and is a 
gapless spin liquid with quasi long range order (algebraically decaying correlations). This behaviour 



persists up to a critical point |37| J2C — 0.241167Ji where a Kosterlitz-Thouless (KT) transition occurs 



to a gapped, dimer phase with exponentially decaying correlations. A simple dimer wavefunction 
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becomes exact at the MG point (J2 = J\/2) f33|- Here it can be shown that the gap is non-zero [|dj. 
At J2 — 00, the lattice separates into two uncoupled Heisenberg chains and is thus again critical. 

For the LS model, the Heisenberg point is not critical, there being a non-zero gap and exponentially 
decaying correlations 0,0- At the Affleck point (J2 = Ji/3) a simple valance bond solid (VBS) 
wavefunction provides the exact ground state f4l|| and the model is rigorously known to be gapped 
Q. The model is solvable at the LS point (J2 = J\) [ p6[ , having soft mode excitations at k — ±27r/3. 
The best numerical studies indicate that this point is critical (i.e. J2 C — Ji), separating a gapped 
VBS (Haldane) phase from a gapless phase with a quasi long ranged 3-fold ground state periodicity 
which persists up to the ferromagnetic boundary ( J 2 = 00), the transition again being of the KT type. 

Until the advent of the DMRG, less was understood regarding the nature of incommensurate 
correlations in these models. In the classical limit S = 00, these models both exhibit a crossover 
to incommensurate behaviour at some threshold J2c\- For the MG model we have J2 C i = Ji/4 and 
^2ci = Ji/2 for the LS model. For J2 < J2c\ correlations are commensurate, spins being antialigned 
in a Neel ground state. For J2 > </2ci the classical pitch angle 6 C \ between successive spins is given by 
C \ = cos _1 ( — Jzcl/Js)- The DMRG has allowed a characterisation of the incommensurate behaviour 
through the accurate calculation of correlations over large distances within large chains. 

In a series of studies of the structure factor 

1 N 

5 ( fc ) ^nT, Ca-e^-i ) fe (14) 
i,j'=i 

where 

C^ = (Vg|5^^g> (15) 

is the correlation function and N the lattice size, a threshold for the onset of incommensurate corre- 
lations J2L was defined as the value of J2 at which the peak in the structure factor A; max begins to 



deviate from ir |10|-|13 . Plots of fc max as well as 6 C \ for the two models are given in Fig. 1. 

The DMRG permits accurate pinpointing of the threshold. For the MG model we have J2L ~ 
0.5206Ji UJ and for the LS model J 2 l ~ 0.4684Ji ffij. Schollwock et al Jill studied C(r), the 
correlation in real space. They observed the existence of a disorder point J2D marking the onset of 
incommensurability in real space. They interpreted J21, as being a Lifshitz point. In fact, it was shown 
for the LS model that J2D = Ji/3 i.e. the disorder point was precisely the Affleck point. Bursill et al 
jL2j showed that for the MG model, J^d — Ji/2 i-e. the disorder point coincides with the MG point. 
J2D lies in the gapped phase, near to where the gap in maximal 111]], |35|| . Numerically it appears as 
though the correlation length £ takes on its minimum at J 2 = J2D- Schollwock et al have pointed out 
that J2ci, J2L and J2D should be expected to merge in the classical limit S — > 00. 



4 The DMRG and critical phenomena 

A critical point in a quantum lattice model is a point in the parameter space at which the ground 
state energy is non-analytic. Generally, the energy gap A (between the ground and first excited states) 
vanishes as the critical point is approached and, correspondingly, the correlation length £ diverges. 
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As a result of the long range correlations in the vicinity of critical points, attempts to characterise 
them using exact diagonalisation methods may be hindered by limitations on the system sizes that 
can be reached. In principle, the DMRG can study very large systems with high accuracy and so it is 
natural to apply the DMRG to the study of critical points. In this section we review attempts to do 
so in magnetic models. We distinguish between methods which characterise critical points via direct 
calculations of A or £, and those which do so by generating renormalisation group transformations in 
the space of coupling constants. 

Direct calculations of A or £ 

The direct calculation of A or £ in the vicinity of a critical point within the DMRG is hindered by the 
fact that the convergence of calculated eigenvalues or correlation functions with m is slowest when A 
is small or £ is large ||, [[ll]. Nevertheless, such calculations have been forthcoming. 

Chitra et al §J| have calculated A for the MG model (|l|) and £ 0, 0, and A have been 
calculated for the LS model ( p"3| ) but in all cases, though high accuracy is possible near the disorder 
point Ji = J2D) results become too inaccurate as the critical point Ji Q is approached for any recovery 
of the finite size scaling results of |37j and Q . 

For the MG model, Bursill et al [|l2| have also defined a dimer order parameter 

D= lim I Cat/2- 1 N/2 - Cjv/2 JV/2+1 1 (16) 

N — >oo 

which characterizes the dimer phase J2 C < J2 < 00, vanishing as J2 — > (^2c) + ■ A plot of D versus 
J2I J\, calculated with m = 200 is given in Fig. 2. Again, the critical point is not pinpointed with this 
approach. 

These somewhat disappointing results are not surprising given that the phase transitions are of the 
KT type, from a gapped phase to an extended gapless phase with an essentially singular divergence 
in £ which is extremely difficult to recover numerically without recourse to scaling theory. 

For transitions such as the Ising transition where the order parameter vanishes only at the critical 
point and the system is gapped elsewhere, there is more hope for direct approaches. The exactly 
solvable 5=1/2 transverse Ising (ITF) model 

h=-jJ2s:-Y, s ? s ! + i ( l7 ) 

i i 

has an Ising transition at j c = 2 which is accurately recovered by the direct approach [42| . 

Kato and Tanaka |l4j have used the DMRG to pinpoint such a transition in the dimerised S = 1 
chain 

H = ^(l-(-#i4i (18) 

i 

at 6 = 0.25 ± 0.01. Also, Pati et al |L5| appear to have discovered such a point in the S = 1 case 
of ( |I2] ) at J2 = 0.730 ± 0.05. Ladder models consisiting of two coupled Heisenberg chains have been 
studied fL6f . The phase transition whereby a gap opens on the introduction of interchain exchange, 
has been characterised. 
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Finally, by implementing new efficiency measures into the finite lattice algorithm, White and 
Affleck |35| have performed calculations on the MG model (|f2|) for lattices with many thousands of 
sites using large values of to. These calulations afford an accurate study of fc max , £ and A for relatively 
large values of J2 . With the help of field theory, this allows a characterisation of the critical point at 
J 2 = 00. 

Coupling constant transformations from the DMRG 

As mentioned, critical points and exponents can be calculated by generating RG transformations in 
the space of coupling constants. This approach has been applied extensively using the RSRG to derive 
the transformations p3[ . 

The idea is that a system of N sites, with Hamiltonian Ti. paramterised by coupling constant (s) 7, 
is divided into N/b blocks of b sites. N/b- fold tensor products of the to lowest energy eigenstates 

B = {\1)i):i = l,...,m} (19) 

of the block Hamiltonians are used to form a truncated basis B ® . . . <g> B (of dimension m N / b ) for the 
system. B is then identified as being a basis for a small repeat unit in a lattice model — a single site, 
or a multisite. For example, if m = 2 then B is identified as the basis for a site containing an S = 1/2 
operator, if to = 4 then identification with a multisite consisting of two S = 1/2 atoms is appropriate. 
Tit ', the projection of TC onto B ® . . . ® B can then be identified with a spin Hamiltonian. 
If H' is of the same form as TL, but with a different constant (vector) 7' then 

7^7' = T(b\j) (20) 
constitutes a RG transformation from which critical points 

7c = T(6| 7c ) (21) 



and critical exponents such as the thermal exponent 

log& 



(22) 



log|T'(6| 7c )| 
can be extracted. 

The RSRG results for the ITF model were quite acceptable [Q. However, for the anisotropic 
5 = 1/2 Heisenberg (XX Z) model 



n 



E [ S * S i+i + \ + h.c.)] (23) 



which is essentially singular [y = 00) at the Heisenberg point (7 C = 1), the RSRG result was incorrect 
[fl5[ in that v \ 2 as b — > 00. 

It is natural to apply the DMRG philosophy to generate coupling constant transformations. That 
is, instead of using the m lowest eigenstates of the block Hamiltonian to form B, the to most probable 
states of a block immersed in an enviroment are used Eg] . This lead to dissapointing results for the 
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ITF model, it being concluded that the DMRG offered no special advantage over the RSRG when 
applied to critical phenomena. 

The studies in jif} can not, however, be considered exhaustive as the block sizes used were quite 
small. In ]47[ , it was proposed that the robustness, portability and systematic improvability of the 
DMRG could be applied to the generation of RG transformations ( po| ) with large blocksizes b, using 
the fact that accurate representations of important states can be gleaned for large blocks. 

This approach was applied to the XX Z model. The result for v : given as a function of b in Fig. 3, 
shows a marked improvement over the RSRG result, the essential singularity being recovered in the 
limit b — + oo p^J . It would seem that more testing of the DMRG as a generator of RG transformations 
might be worthwhile. 

Finally, it was noted in p7| that the values of the important eigenvalues in the various symmetry 
sectors of the density matrix provide a general and intuitive scheme for choosing a suitable value of m 
and hence for identifying the renormalised Hamiltonian Ji' . For instance, it was noted that the 5* = 1 
Heisenberg model mapped naturally onto a ladder model under this scheme, as has been demonstrated 
analytically and numerically p§[ ]. 

5 Recent extensions 

In this section we discuss some current trends in advancing the applicability of the DMRG. We focus 
on recent extensions of the DMRG to the calculation of thermodynamical properties of 2D classical 
systems and (hence) ID quantum systems, as well as to the calculation of low energy properties of 
quantum models in 2D. 

2D classical and ID quantum systems at non-zero temperature 

In a recent paper, Nishino exploited the similarities between the transfer matrix of a 2D classical 
system and the Hamiltonian of a ID quantum system in order to formulate the DMRG for 2D classical 
systems ^9) . Results for the specific heat of the 2D Ising model are most promising, the exact result 
being recovered with good accuracy over the whole temperature range, even at criticality, using a 
reasonably modest basis set jl9| . 

Potential applications of this approach include the determination of the critical properties of 2D 
classical models with high spin, anisotropy and frustration. Also, this extension paves the way for the 
study of quantum spin chains at non-zero temperature, as the Trotter-Suzuki decomposition allows 
such systems to be mapped onto classical 2D systems. 

This programme has been carried out for the dimerised XY chain 

n = - L S* + Sl^) - 7 £(^S 2 * i+1 + SlSl +1 ) (24) 

i i 

Although technical difficulties with the transfer matrix limit the basis sizes which can be used, very 
reasonable results, down to low temperatures T have been obtained for the internal energy u, especially 



for 7 7^ 1) when the model possesses a substantial gap [ p0| . A plot of u versus T, for 7 = 2 is given 
in Fig. 4. 

Finally, Nishino and Okunishi have derived two further reformulations of the DMRG — the product 
wavefunction renormalisation group (PWFRG) |5lj] and the corner transfer matrix renormalisation 
group (CTMRG) p^l . These formulations may offer dramatically improved efficiency over the DMRG 
and the means of calculating dynamical correlation functions in spin chains |32]]. The computational 
efficiency of the CTMRG has been exploited so as to obtain highly accurate results for the 2D Ising 
model at criticality |52]| . 

Two dimensional quantum systems 

Soon after the inception of the DMRG, attempts were made at generalising the finite lattice algorithm 
to two dimensions |5^]-p^|. However, various conditions which permit the unprecedented accuracy of 
the DMRG in ID — the fact that, for short ranged interactions there is no interaction between system 
and environment blocks, as well as the preservation of reflection symmetry, do not persist in higher 
dimensions. This is reflected in the relatively slow convergence of the DMRG when applied without 
modification to two dimensional systems | p3[ . 

In order to circumvent this problem, the DMRG has been reformulated for fermion systems in 
terms of momentum space basis states p5| . This advancement is made possible through an elegant 
solution to the technical problem of storing an updating the potential energy operator which is a sum 
of products of four creation and annihilation operators. The formulation offers the following distinct 
advantages: 

1. The total momentum can be used as a quantum number in diagonalising the supcrblock Hamil- 
tonian. 

2. Rotational symmetry is not broken. The breaking of this symmetry is a major cause of the slow 
convergence of the real space formulation of the DMRG. 

3. A single program structure applies quite generally in all dimensions. 

Results for the 2D Hubbard model are extremely promising |55| . Agreement with exact calculations 
on 4 x 4 lattices is good as is agreement with projected quantum Monte Carlo (PQMC) and stochastic 
diagonalisation (SQ) data for 8x8 lattices. In fact, the DMRG, being a truncated basis expansion, 
generates a variational bound for the ground state energy. This bound is better than that obtained 
from SQ and is thus the best variational bound on the ground state energy to date (see table 1). Also, 
a numerical solution of the Hubbard model on a 12 x 12 lattice has been afforded by this scheme [ p5| . 

Potential applications of this approach include heavily doped fermion systems, fermion systems 
with weak coupling and fermion systems with long range coulomb interactions. 

Unfortunately, spin hamiltonians or models with dressed fermion operators such as the t-J model 
cannot be conveniently expressed in terms of momentum space operators. However, there have been 
algorithmic advancements in the DMRG which have been exploited to permit the study of such 
models in 2D |37j-[|5£j. White has performed successful calculations of the spin gap in a (frustrated) 
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N 


DMRG 


SQ a 


10 


-34.325 


-34.31 


18 


-54.394 


-54.37 


26 


-66.098 


-66.05 




a) Ref. E 





Table 1: Comparison of DMRG and SQ results for the ground state energy of the 2D Hubbard model 
on an 8 x 8 lattice, with 11 = 41 in the usual notation, for various electron fillings N. 

2D Heisenberg model for the spin-gapped system CaV40g using a 24 x 11 lattice J57| and White and 
Scalapino have investigated the nature of one and two hole ground states in the t-J model on 10 x 7 
clusters |^8| . 

These calculations have been made possible by the use of highly efficient code which minimizes the 
number of hamiltonian matrix multiplications required in the sparse diagonalisation. Moreover, the 
problem of symmetry breaking which leads to the need for a large number of matrix multiplications 
has been recently addressed J59| . 

In |5S|1 an optimal scheme for sweeping the lattice has been developed. This has afforded the study 
of a frustrated Heisenberg model (i.e. with next-neighbour exchange) on a 20 x 20 lattice. Good 
convergence along with the accurate recovery of exact results on the 4x4 lattice permits confident 
characterizations of the phases for various values of the coupling constants. 
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Figure cations 

1. DMRG results for the structure factor peak position A: max and classical pitch angle 6 C \ (in units 
of 7r) as functions of the frustration J\/ J2 for two frustrated spin systems — the MG and LS 
models. 

2. DMRG results for the dimcr order parameter D as a function of the frustration J1/J2 for the 
MG model. The known position of the phase transition is indicated by the dashed line. 

3. DMRG result for the critical exponent v as a function of blocksize b for the anisotropic Heisenberg 
model. The dashed line is the result from the RSRG approach. 

4. DMRG result (dashed line) for the internal energy as a function of temperature for the dimerised 
XY model with dimerisation ratio 7 = 2. The full line is the exact result. 



13 



